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Abstract. Shearlet tight frames have been extensively studied during the last years due 
to their optimal approximation properties of cartoon-like images and their unified treatment 
of the continuum and digital setting. However, these studies only concerned shearlet tight 
frames generated by a band- limited shearlet, whereas for practical purposes compact support 
in spatial domain is crucial. 

In this paper, we focus on cone-adapted shearlet systems which - accounting for stabil- 
ity questions - are associated with a general irregular set of parameters. We first derive 
sufficient conditions for such cone-adapted irregular shearlet systems to form a frame and 
provide explicit estimates for their frame bounds. Secondly, exploring these results and 
using specifically designed wavelet scaling functions and filters, we construct a family of 
cone-adapted shearlet frames consisting of compactly supported shcarlets. For this family, 
we derive estimates for the ratio of their frame bounds and prove that they provide optimally 
sparse approximations of cartoon-like images. 



1. Introduction 

Over the last 20 years wavelets have established themselves as a key methodology for 
efficiently representing signals or operators with applications ranging from more theoretical 
tasks such as adaptive schemes for solving elliptic partial differential equations to more 
practically tasks such as data compression. It is fair to say that nowadays wavelet theory 
can be regarded as an essential area in applied mathematics. A major breakthrough in this 
field was achieved through Daubechies in 1988 by introducing compactly supported wavelet 
orthonormal bases with favorable vanishing moment and support properties [3]. 

Recently, it was proven by Candes and Donoho in pTJ that wavelets do not perform op- 
timally when representing and analyzing anisotropic features in multivariate data such as 
singularities concentrated on lower dimensional embedded manifolds, where their focus was 
on the 2-dimensional situation. This observation initiated a flurry of activity to design a 
representation system which provides optimally sparse approximations of so-called cartoon- 
like images while having as favorable theoretical and computational properties as wavelet 
systems (see, e.g., [HE])- In 2005, shearlets were developed by Labate, Weiss, and two of 
the authors in [2U] (see also [12]) as the first directional representation system with allows 
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a unified treatment of the continuum and digital world similar to wavelets, while providing 
(almost) optimally sparse approximations within a cartoon-like model [33] - 'almost' in the 
sense of an additional log-factor which is customarily regarded as negligible. 

Several constructions of discrete shearlet frames are already known to date, see [T2j [161 
[21 [2U [15]. However, taking applications into account, spatial localization of the analyzing 
elements of the encoding system is of uttermost importance both for a precise detection of 
geometric features as well as for a fast decomposition algorithm. But, one must admit that 
non of the previous approaches encompasses this crucial case. Hence the main goal of this 
paper is to provide a comprehensive analysis of cone-adapted - the variant of shearlet systems 
feasible for applications - discrete shearlet frames encompassing in particular compactly 
supported shearlet generators. Our contribution is two-fold: We firstly provide sufficient 
conditions for a cone-adapted irregular shearlet system to form a frame for L 2 (R 2 ) with 
explicit estimates for the ratio of the associated frame bounds; and secondly, based on these 
results, we introduce a class of cone-adapted compactly supported shearlet frames, which 
are even shown to provide (almost) optimally sparse approximations of cartoon-like images, 
alongside with estimates for the ratio of the associated frame bounds. 

1.1. Shearlet Systems. Referring to the detailed introduction of shearlet systems in Sec- 
tion [21 we allow us here to just briefly mention those main ideas and results, which are crucial 
for an intuitive understanding. 

Shearlet systems are designed to efficiently encode anisotropic features. In order to achieve 
optimal sparsity, shearlets are scaled according to a parabolic scaling law, thereby exhibiting 
a spatial footprint of size 2~i times 2~- J / 2 . They parameterize directions by slope encoded 
in a shear matrix, since choosing shear rather than rotation is in fact decisive for a unified 
treatment of the continuum and digital setting. We refer for more details to [121 EZ] for the 
continuum theory and [191 HQl 121] for the digital theory. 

Shearlet systems come in two ways: One class being generated by a unitary representation 
of the shearlet group and equipped with a particularly 'nice' mathematical structure, how- 
ever causes a biasedness towards one axis, which makes it unattractive for applications; the 
other class being generated by a quite similar procedure however restricted to a horizontal 
and vertical cone in frequency domain, thereby ensuring an equal treatment of all direc- 
tions. For both cases, the continuous shearlet systems are associated with a 4-dimensional 
parameter space consisting of a scale parameter measuring the resolution, a shear param- 
eter measuring the orientation, and a translation parameter measuring the position of the 
shearlet. A sampling of this parameter space leads to discrete shearlet systems, and it is 
obvious that the possibilities for this are numerous. A canonical sampling approach - using 
dyadic sampling - leads to so-called regular shearlet systems. However, lately, questions 
from sampling theory and stability issues brought the study of irregular shearlet systems to 
the researcher's attention. 

1.2. Necessity of Compactly Supported Shearlets. Shearlet theory has already im- 
pacted various applications for which the sparse encoding or analysis of anisotropic features 
is crucial, and we refer to results on denoising [21], edge detection [2], or geometric sep- 
aration [H [9]. However, similar as in wavelet analysis, a significant improvement of the 
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applicability of shearlets is possible, if compactly supported shearlet frames would be made 
available. Let us discuss two examples to illustrate the necessity of such a study. 

Imaging Science. In computer vision, edges were detected as those features governing an 
image while separating smooth regions in between. Thus imaging tasks typically concern 
an especially careful handling of edges, for instance, avoiding to smoothen them during a 
denoising process. Therefore, when exploiting a decomposition in terms of a representation 
system for such tasks, a superior localization in spatial domain of the analyzing elements 
is in need. This would then allow a very precise focus on edges, thereby reducing or even 
avoiding various artifacts in the decomposition and analysis process. Secondly, due to the 
increased necessity to process very large data sets, fast decomposition algorithms are crucial. 
Learning from the experience with structured transforms such as the fast wavelet transform, 
it is conceivable that a representation system consisting of compactly supported elements 
will provide a significant gain in complexity. 

Partial Differential Equations. Adaptive schemes for solving elliptic partial differential 
equations using wavelet decompositions have turned out to be highly beneficial both in 
theory and applications (cf. [3]). However, hyperbolic partial differential equations exhibit 
shock fronts, which are not optimally sparsely encoded by wavelets due to their isotropic 
footprints. For a directional representation system such as shearlets to be feasible as an 
encoding method for an adaptive scheme, one main obstacle is the sparsity of the associated 
stiffness matrix. However, with good localization properties such as compact support, this 
obstacle might be overcome. Another problem is the handling of boundary conditions on 
a bounded domain, which seems also attackable by having compactly supported shearlet 
frames at hand. Last, but not least, again a fast decomposition algorithm is of importance, 
and we refer to the previously described application for a reasoning why having compact 
support will make a tremendous difference here. 

1.3. Previous Work on Shearlet Constructions. Let us now discuss the history and 
the state-of-the-art of shearlet constructions. The first class of shearlets, which were shown 
to generate a tight frame, were band-limited with a wedge-like support in frequency domain 
specifically adapted to the shearing operation (see [201 12] )• This particular class of cone- 
adapted shearlet frames was already extensively explored, for instance, for analyzing sparse 
approximation properties of the associated shearlet frames [T3] . 

Shortly afterwards, a different avenue was undertaken in [16J, where a first attempt was 
made to derive sufficient conditions for the existence of irregular shearlet frames, which, 
however, could not decouple from the previous focus on band-limited shearlet generators. In 
addition, this result was stated for shearlet frames which arose directly from a group repre- 
sentation. In some sense, this path was continued in [2], where again sufficient conditions for 
this class of irregular shearlet frames were studied. Let us also mention at this point that 
a more extensive study of these systems was performed in [15] with a focus on necessary 
conditions and a geometric analysis of the irregular parameter set. 

Now returning to the situation of cone-adapted shearlet frames, a particular interesting 
study was recently done in [2T], where a fast decomposition algorithm for separable shearlet 
frames - which are non-tight - was introduced, however without analyzing frame or sparsity 
properties. We should mention that it was this work which gave us the intuition of taking 
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a viewpoint of separability for the construction of compactly supported shearlet frames in 
this paper. The question of optimal sparsity of a large class of cone-adapted compactly 
supported shearlet frames was very recently solved by two of the authors in [18]. However, 
no sufficient conditions for or construction of cone-adapted compactly supported shearlet 
frames with estimates for the frame bounds nor with sparsity analyses are known to date. 
This is an open question even for other directional representation systems such as the perhaps 
previously most well-known curvelet system pp. 

1.4. Our Contribution. Our contribution in this paper is two-fold. Firstly, Theorem 13.41 
provides sufficient conditions on the irregular set of parameters and the generating shearlet 
to form a cone-adapted shearlet frame alongside with estimates for the ratio of the frame 
bounds. This result greatly extends the result in [16] by, in particular, encompassing com- 
pactly supported shearlet generators. Secondly, exploring these results and using specifically 
designed wavelet scaling functions and filters, we explicitly construct a family of cone-adapted 
shearlet frames consisting of compactly supported shearlets with estimates for their frame 
bounds in two steps: For functions supported on the horizontal cone in frequency domain 
and then for functions in L 2 (M. 2 ) (see Theorems 14.71 and 14.91) . This construction is the first of 
its kind in shearlet theory, and it is worth mentioning that, to the author's knowledge, even 
in wavelet theory there does not exist a comparable one for compactly supported wavelets. 
Theorem 14.101 then proves that this family does even provide (almost) optimally sparse ap- 
proximations of cartoon-like images. 

1.5. Outline. In Section &\ we introduce the necessary definitions and notation for cone- 
adapted irregular shearlet systems. In particular, in Subsections 12.11 and 12.21 we present the 
viewpoint of deriving cone-adapted irregular shearlet systems from sampling the parame- 
ters of cone-adapted continuous shearlet systems for the first time. Subsection 12.31 is then 
concerned with introducing a large class of cone-adapted irregular shearlet systems, coined 
feasible shearlet systems, which will be the focus of this paper. Sufficient conditions for such 
feasible shearlet systems to form a frame and explicit estimates for their frame bounds are 
then discussed in Section |3j In Section HI these results are explored to construct a family 
of cone-adapted regular compactly supported shearlet frames with estimates for the ratio of 
the associated frame bounds (Subsections 14. II and 14.21) . In Subsection 14.31 h is finally proven 
that this family does even provide (almost) optimally sparse approximations of cartoon-like 
images. Some very technical, lengthy proofs are deferred to Section 

2. Regular and Irregular Shearlet Systems 

2.1. Cone-Adapted Shearlets. Shearlets are highly anisotropic representation systems, 
which optimally sparsify C 2 (lR 2 )-functions apart from C 2 -discontinuity curves. Their main 
advantage over other proposed directional representation systems is the fact that they provide 
a unified treatment of the continuum and digital world. Today, shearlets are utilized for 
various applications, and we refer to [TU El E] as some examples. 

Introduced in [TJ], discrete shearlet systems - shearlet systems for L 2 (R 2 ) with discrete 
parameters - exist in two different variants: One coming directly from a group representation 
of a particular semi-direct product, the so-called shearlet group, and the other being adapted 
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to a cone- like partitioning of the frequency domain. In fact, the second type of shearlet 
systems are the ones exhibiting the favorable property of treating the continuum and digital 
setting uniformly similar to wavelets, and are the ones relevant for applications. Hence those 
are the discrete shearlet systems we focus on in this paper. 

We now introduce cone-adapted discrete regular shearlet systems - previously also referred 
to as 'shearlets on the cone' - as a special sampling of the (cone-adapted) continuous shearlet 
systems, which were exploited in [IT] . Since all previous papers consider discrete regular 
shearlet systems as sampling the parameters of continuous shearlet systems refer to the first 
type of shearlet systems arising from a group representation (see, e.g., [161 12]), this is the 
first time this view is taken and hence is presented in a more elaborate manner. 

2.1.1. (Cone- Adapted) Continuous Shearlet Systems. Shearlets are scaled according to a par- 
abolic scaling law encoded in the parabolic scaling matrices A a or A a , a > 0, and exhibit 
directionality by parameterizing slope encoded in the shear matrices S s , s G M, defined by 

^=(o Va) OT A "=(\f a) 

and 

Ss = (o l) ' 

respectively. 

To ensure an (almost) equal treatment of the different slopes, which is evidently of ut- 
termost importance for practical applications, we partition the frequency plane into the 
following four cones C\ - C^. 

[ {(£i,6)eM 2 :£i>l, 16/61 <1} : L = h 
r= \ {(6,6)eM 2 :6>l, 16/61 <1} : ' = 2, 

L I {(6,6) e M 2 :6<-1, 16/61 < 1} = t = 3, 
{ {(6,6) e R 2 :&<-!, 16/61 < 1} = t = 4, 

and a centered rectangle 

7i = {(6,6)eM 2 :||(6,6)l|oo<i}. 

For an illustration, we refer to Figure Ufa). 

The rectangle 1Z corresponds to the low frequency content of a signal, which is customarily 
represented by translations of some scaling function. Anisotropy now comes into play when 
encoding the high frequency content of a signal, which corresponds to the cones C\ - C4, 
where the cones C\ and C3 as well as C2 and C4 are treated separately as can be seen in the 
following definition. 

Definition 2.1. The (cone-adapted) continuous shearlet system SW((j),ip,ip) generated by 
a scaling function <fi 6 L 2 (M. 2 ) and shearlets 6 L 2 (M?) is defined by 

su(<t>, = $(0) u$(f)u vj>o/o, 

where 

$(0) = {0 t = 0(--t):tGM 2 }, 
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(a) (b) 

Figure 1. (a) The cones C\ - C 4 and the centered rectangle 1Z in frequency 
domain, (b) The tiling of the frequency domain induced by discrete shearlets 
with the support of one shearlet from Example 12.31 exemplarily highlighted. 

*(V0 = &a,s,t = a-i^A-'S^i - - 1)) : a G (0, 1], s G [-(1 + a 1 / 2 ), 1 + a 1 / 2 ], t G I 
and 

= {4>a,s,t = a"f ^(i a " 1 5 s - T ( ■ -t)):ae (0, 1], s G [-(1 + a 1 / 2 ), 1 + a 1 / 2 ], t G I 
where a indexes scale, s indexes shear, and t indexes translation. Setting 

Scone = {(a,s,t) :aG (0,1], sG [-(1 + a 1/2 ), 1 + a 1 ' 2 }, t G M 2 }, 



the associated (Cone-Adapted) Continuous Shearlet Transform S'K^^^f 



of some function f G L 



is given by 



Notice that in [T7] only special generators ^ were considered, whereas here we state the 
definition for all ip G L 2 (M 2 ); see also the interesting extension [TT] . 



2.1.2. f Cone-Adapted) Discrete Shearlet Systems. A discretization of the shearlet parameters 
is typically achieved through sampling of the parameter set of (cone-adapted) continuous 
shearlet systems, where 1R 2 is sampled as c{E? and, for the horizontal cone C\ U C3, E> cone is 
discretized as 

{(2-'', k2-il\S k2 - ]/ ,A 2 - ] M c m : 3 >0,ke {-\7?'% 

where, for c = (ci, c 2 ) G 



\2^ 2 ]}, mGZ 2 }, 



i + ) 2 , M c denotes the sampling matrix 

M„ 



ci 
.0 c 2y 

For the vertical cone C 2 U C 4 , § co „ e is discretized in a similar way, now using the sampling 
matrix 

'Ca 0' 

Cl/ 

The anisotropic sampling of the position parameter is chosen to allow the flexibility to 
oversample only in the shearing direction, which will reduce the redundancy of the generated 
system. This will later become apparent in the ratio of the associated frame bounds. 
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Also the range of values for k deserves a comment, since the two values — 2 J//2 — 1 and 
2 j / 2 + 1 which arose from the sampling procedure were replaced by - \2^ 2 ] and \2^ 2 ] . The 
reason for this is that 2 J / 2 is not an integer for odd scales j. Also, for the classical almost- 
separable shearlet generator introduced in Example 12.31 the support of the Fourier transform 
of the associated shearlets does not anymore fall into CiUC 3 for k = -2^ 2 -l and k = 2^ 2 + l 
(if considering the horizontal cone). 

This now gives the following discrete system. 

Definition 2.2. For some sampling vector c = (01,02) G (R + ) 2 , the (cone-adapted) regu- 
lar discrete shearlet system S7i(c;(f),xjj,ip) generated by a scaling function <p G L 2 (R 2 ) and 
shearlets G L 2 (M 2 ) is defined by 

SH(c; 0, V, = $(ci, <P) U tt(c, V) U §(c, $), 

where 

$(ci, <p) = {<p rn = <j>(- - c x m) : m G Z 2 }, 
*(c,V) = {^,fc,m = 2 3 ^V(,S_ fc A 2J ■ -M c m) : j > 0, < [2^ 2 ],m G Z 2 }, 

and 

tt(c,$ = {^ ijfcjro = 2 Zj '^(S T _ k A 23 ■ -M c m) : j > 0, < [2^ 2 ],m G Z 2 }. 

Setting 

Acone = {(j, k, m):j>0,\k\< [2 j/2 ] , m G Z 2 }, 
i/ie associated (Cone- Adapted) Regular Discrete Shearlet Transform SW^^^f : Z 2 x A 2 one — > 
C 3 of some function f G L 2 (M 2 ) is given fry 

SK^^fi 171 '' 0'> m )' 07 ™)) = ((/» 0m' )> (/» fern), (/, 

The reader should keep in mind that although not indicated by the notation, the functions 
(frm, il>j,k,mi an d ipj,k s m & U depend on the sampling constants c\,c<i- Notice further that the 
particular sampling set (JT]) forces a change in the ordering of parabolic scaling and shearing, 
which will be mimicked later by the class of irregular parameters we will be considering. The 
definition itself shows that the shearlets ipj,k,m (also ipj^,m) live on anisotropic regions of size 
2 -J ' x 2~2 at various orientations in spatial domain and hence on anisotropic regions of size 
2 J x 2a in frequency domain (see also Figure DJb)). 

Example 2.3. The classical example of a generating shearlet is a function ip G L 2 (K 2 ) 
satisfying 

^(0 = ^(6,6) = Mi) 4(g), 

where ipi G L 2 (M) is a discrete wavelet, i.e., satisfies the discrete Calderon condition given 
by E jeZ l4(2-y)| 2 = 1 for a.e. cue R, with & G C°°(R) and supp & C [-|, -±] U [i, §], 
an<i ^2 G L 2 (M) is a 'bump' function, namely Y2l=-i + ^)| 2 = 1 / or a - e - w6 [-1,1], 
satisfying tp 2 G C°°(R) and supp^ 2 C [—1,1]. There are several choices of and ip2 
satisfying those conditions, and we refer to [12] for further details. 

The tiling of frequency domain given by this band-limited generator and choosing ip = 
ijj(x2, X\) is illustrated in Figure\^b) . For this particular choice, using an appropriate scaling 
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function forlZ, it was proven in [121 Thm. 3] that the associated (cone- adapted) discrete 
regular shearlet system S%{1; 0, ip, ip) forms a Parseval frame for L 2 (IR 2 ). 

2.2. Irregular Sampling of Shearlet Parameters. Questions arising from sampling the- 
ory and the inevitability of perturbations yield the necessity to study shearlet systems with 
arbitrary discrete sets of parameters. To derive a better understanding of this problem, we 
take the viewpoint of regarding the previously considered set of parameters 

{(2~ j : k2~3l\ S k2 - 3/ 2A 2 - 3 M c m : j > 0, k G {-[2^], . . . , [2^]}, m G Z 2 } 

as a discrete subset of S cone (for C\ U C3 and similarly for C\ U C3 with M c substituted by 
M c ) and for the low frequency part M. 2 likewise. Then, an arbitrary set of parameters is 
merely a different sampling set, and questions of density conditions of such a sampling set 
or its geometric properties are lurking in the background. Some density-like properties will 
indeed come into play in this paper mixed with decay conditions on the shearlet generator. 
For an extensive analysis of the geometric properties of the sampling set yielding necessary 
conditions for shearlet frames directly arising from a group representation, we refer to [15]. 
Next we formally define (cone-adapted) irregular shearlet systems. 

Definition 2.4. Let A and A, A be discrete subsets of M? and §> cone , respectively, and let 
G L 2 (M 2 ) as well as ip.^p G L 2 (M 2 ). Then the (cone-adapted) irregular discrete shearlet 
system ST-L(A, A, A; 0, 1/}, ijj) is defined by 

SH(A, A, A; 0, ^ = $(A, 0) U *(A, j>) U *(A, i>), 

where 

$(A,0) = {0 t = 0(--t):tGA}, 
*(A,V) = {i>a, s ,t = a-i^A-'S^i - -t)) : (a,s,t) G A}, 

and 

(A, $) = {4> a , Stt = a--4(A- 1 S; T ( ■ - t)) : (a, s, t) G A}. 
Then the associated (Cone- Adapted) Irregular Discrete Shearlet Transform SH^^^f : A x 
A x A — y C 3 of some function f G L 2 (IR 2 ) is given by 

2.3. Class of Irregular Shearlet Systems. Since the low frequency part is already very 
well studied, and since both cone pairs C\ and C3 as well as C 2 and C4 are treated similarly, 
we from now on restrict our attention to the horizontal cone 

C :=dUC 3 , 

and are interested in frame properties of the system \&(A, ip), when expanding functions in 

L 2 (C) = {/GL 2 (M 2 ):supp/CC}. 

We remark that it is sufficient to consider the function space L 2 (C), since in classical shearlet 
theory, the representing shearlets are orthogonally projected onto the cones (and, in partic- 
ular, onto the cone C under consideration). For clarity, we further introduce the notion 

SH(A,i>) :=¥(A,^). 
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We will now focus on a class of sets of parameters satisfying some weak conditions, which 
in a sense endow them with - loosely speaking - similar properties as ([T]). We remark that, 
since one main motivation for considering irregular sets of parameters are stability questions, 
such a constraint seems very natural. Also, we impose weak decay conditions on the shearlet 
generator, which still include spatially compactly supported shearlets; our main objective. 

2.3.1. Feasible Sets of Parameters. We now first discuss our hypotheses on the set of param- 
eters. For this, let 

{{aj,s j)k ,f jtkt J :j>0,kE Kj,m E Z 2 }, K 3 CZ,c = ( Cl ,c 2 ) E (M + ) 2 , 

be an arbitrary discrete set of parameters in § cone . 

For the decreasing sequence {aj}j>o C M + to act as a scaling parameter, we require a 
certain growth condition, which we phrase as follows: For each \l E (0, 1), there shall exist 
some positive integer p such that 

^ < /a for all j > 0. (2) 

a j 

Obviously, this condition is satisfied for the customarily utilized sequence aj = with 
p = 1 and fi = 1/2. 

Further, we assume the sequence {sj s k}j,k C K of shear parameters to be of the form 

Sj,k = Sky/a] for each j > 0, k E Kj, 

which allows a change of the parabolic scaling and shear operators similar to the discretization 
of the (cone-adapted) continuous shearlet systems (see Definition 12. 2j) . We also restrict to 

|sfe| < a 1 / 2 + 1 for all j > 0, k E Kj, (3) 

and mimic the discretization in the regular case. From now on, we assume that E Kj for 
each j > and Sq = for simplicity, and remark that all of the results in this paper can 
be easily extended to a general sequence of shear parameters Sk without this assumption. 
Finally, we assume the sequence {sk}k to satisfy a density-like condition, which will ensure 
that the shearing operation does not cause a 'piling up' of the essential supports of the 
shearlets. In particular, we require that 

sup ^min{l, \(S Sk x) 2 \} ■ min{l, \(S Sk x)i\~' y } < C(j) < oo for any 7 > 1. (4) 

x=(xi,X2)eM 2 fcgZ 

The sequence {t C j :k m }j t k,m C M 2 of translation parameters shall have a grid structure, which 
we impose by assuming that 

tj,k,m = S Sj k A aj M c m, for some c = (ci, c 2 ) E (M + ) 2 . (5) 
Summarizing, we get the following 

Definition 2.5. Let the sequences {Oj}j>o C M + ; { s j,k}j>o,keKj C K, and the sequence 
{Pj,k m}j>0MKj,m€Z 2 C M 2 satisfy (|2j) ; (J3j) — (SJ) ^ and (jSJ), respectively. Then we call the se- 
quences 

{(aj, s jtk = s ky /a~j,t C j k m = S Sj k A a M c m) : j >0,k E Kj,m E Z 2 } 
a feasible set of parameters. 
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2.3.2. Feasibility Conditions for Generating Shearlet. We next impose the following decay 
condition on the generating shearlets ip £ L 2 (M 2 ): 

Definition 2.6. A function ip £ L 2 (R 2 ) is called a feasible shearlet, if there exist a > 7 > 3, 
and q>q'>0,q>r>0 such that 

h?(£i,6)| < min{l, \q^\ a } ■ min{l, • min{l, H 2 \^}. (6) 

At first sight it seems picky to choose three different constants q, q', and r for the dif- 
ferent 'decay areas' in frequency domain. However, this flexibility will become crucial for 
deriving estimates for the lower frame bound, and hence for analyzing the class of compactly 
supported shearlet frames we aim to introduce. 

The class of feasible shearlets certainly includes the classical example of generating shear- 
lets (Example I2.3j) as well as all L 2 -functions whose Fourier transform is of polynomial 
decay of a certain degree. In particular, it includes functions which are spatially compactly 
supported. Moreover, condition ([6]) is chosen in such a way that ip exhibits an essential 
support in frequency domain somehow similar to the support of the classical shearlets de- 
fined in Example 12.3} but is adapted to separable shearlet frames which later on will serve 
as the structure for our explicitly constructed compactly supported shearlet frames. For an 
illustration we refer to Figure [2J 




Figure 2. Essential support of in the frequency domain for a shearlet ijj 
satisfying the decay condition d^J). 



2.3.3. Class of Shearlet Systems. Concluding, in this paper we are concerned with a partic- 
ular class of (cone-adapted) irregular shearlet systems defined as follows. 

Definition 2.7. Let A = {{o>j,Sj i k,tj krr ^) : j > 0, k £ Kj,m £ Z 2 } be a feasible set of 
parameters and let ip £ L 2 (M. 2 ) be a feasible shearlet. Then we call 

SH{A,*P) = {^ k , m = a~^(S_ Sk A a -i ■ -M c m) :j>0,ke Kj,m £ Z 2 } 
a feasible shearlet system. 
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Notice that such a system depends in particular on the parameters /i, a, 7, q, q', and r as 
well as on the function C(-). We further remark the abuse of notation for if)j t k,m a s compared 
to Definition 12.21 to keep the notation simple. However, it will always be made clear which 
interpretation is meant. 

3. A General Sufficient Condition for Shearlet Frames 

In this section, we will state a very general sufficient condition for the existence of cone- 
adapted irregular shearlet frames. A similar condition can be found in [16], however for 
shearlet frames arising from a group representation. This result can not be directly carried 

—1/2 

over, since the relation between the cone, the bound a- for the parameters s&, and the 
most likely not band-limited shearlet generator if) needs careful handling. Also, the result in 
[T6] does not include compactly supported shearlet generators. 



3.1. Covering Properties. For stating and proving the result, we require a particular 
function, which is an adapted form already appeared in [16] (see also [5] for the wavelet 
situation). It can be regarded as a variant of the main term in the t q condition for tight 
affine frames [5]. 

Let 5"H(A, if)) be a feasible shearlet system as defined in Definition 12. 71 The main objective 
will now be the function $ : C x R 2 — > E defined by 

®(M = EE msJ h A aj mif;(SlA a ^ + u))\. (7) 
j>o keKj 

Notice that this function measures the extent to which the essential supports of the scaled 
and sheared versions of the shearlet generator overlaps. For later use, we also introduce the 
function r : M 2 — > R defined by 

r(a;) = ess sup to), 

measuring the maximal extent to which these versions overlap for a fixed distance, as well 
as the values 

Linf = essinf $(£, 0) and L sup = ess sup 0), (8) 

« eC tec 

which relate to the classical discrete Calderon condition. Finally, we also require the values 
R(c)= i T ( M c lm ) r (-M" 1 ™)] 172 , where c = (c u c 2 ) G (R + ) 2 . (9) 

m£Z 2 \{0} 

The function R(c) measures the average of the symmetrized function values T(M~ 1 m). 

We now first turn our attention to the terms L sup and R{c) and provide explicit estimates 
for those, which will later be used for estimates for frame bounds associated to a shearlet 
system. We start by estimating L sup in 

Proposition 3.1. Let Sl-L(h,il)) be a feasible shearlet system, and let L sup be defined as in 
©. Then 

L, v < p • (2 ■ 0(2,)) ( [ l„g 1/(( r -L rT + r -L 57 ) < 00. (10) 
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In the special case Sk = k, k £ Z, the following stronger estimate holds: 

L„ p < P ■ (I (2 + ^) + l) (f log,,,. (I)] + T -^_ r + j-^) < 00. (11) 

One essential ingredient in the proof of this lemma is a result from [22], which we state 
here for the convenience of the reader. 

Lemma 3.2. [221 Lem. 2.1] Let \i £ (0,1), p £ Z + , and Ze£ {aj}jez C K + 6e a decreasing 
sequence which satisfies (0|). Then, for all a > and £ > 0, 

y a -< p r-^, y fl «<p*«_!_. 

^ J 1 - u Q ^ J ~ 1 - u a 

This now allows us to prove Proposition 13.11 
Proof of Proposition\3Ji Let £ = (6,6) T e M 2 . By (jSJ), we obtain 

$(e,o) 

< sup Vmin{l, |g % £ 1 | 2Q }min{l, IgVil"' 7 } E min {l, |r(aj /2 6 + s fc a 3 -6)|- 2 T}.(12) 
Letting 771 = g£x, 

$(£, 0) < sup ^^min {1, |aj77i| 2a_1 } min{l, |g'g _1 aj?7i| 27 } - min {l, |rg~ 1 aj77i|} 

l/2 ( 



• min{l, |ray £ 2 + rg 27 }- (13) 

By©, 

y^-min{l, |rg~ 1 aj?7i|}min{l, |ra| /2 £ 2 + rq^a^xSkf 21 } < -C(2 7 ), 



k£Z 

Hence we can continue ( {TBI) by 



$(£,0) < f^C(2 7 )) sup Vmin{l,| ai r/ 1 | 2Q - 1 }min{l,| g V 1 a^i|- 27 } 
-C(2 7 )J sup(^X[o ) i)(l a ^il)l a i^i| 2a " 1 + X[i,4)(l a i r Zil) 



i>o 

r '' iCJ1 " j>0 

+X[^,oo)(|aj^i|)|gV 1 «j"ir 27 ^ 

< (^c(2 7 )) ■ sup ( Yl i^ii 2a_1 +Ex[i,*](i^ii)+ E i^v^r 27 

m |oyr/i|<l j>0 | g 'q-ia 3 r7i|>l 

The first claim ( TTOj) now follows from Lemma 13.21 
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Next assume that = k for all fcGZ. To derive an estimate for L sup in this case, we first 
split (fT2~j) into the cases k = and k ^ with r\\ = g£i and 772 = ra 'j^2, which yields 

$(e,o) 

< sup min {1, 



a,?7i| 



j>o 



}min{l,|g'g 1 a J ?7 1 1 27 } sup - min{l, |rg 1 a :7 -?7i | } 

r?26[0,|r-g- 1 a J -77i|] fe ^ Q r 

min{l, 1 772 + rq~ 1 a.jr]ik\~' 2 ' 1 } + sup y^ min {1, |aj?7i| 2a } min{l, \q'q~ 1 aji] 1 \~ 2j }. (14) 



j>0 



Notice that 



sup - min{l, \rq 1 a 7 r? 1 1 } minil, \n 2 + rq 1 a j r]ik\ 27 } < - ( 2 H J 

[o,N-^i|]^ r ~ r V 2 7 -i; 



and 



min {1, |aj?7i| 2a } min{l, \q'q 1 ajfj 1 \ 27 } < min {1, \ajTji\ 2a 1 }min{l,|g'g 1 a J ?7 1 1 27 }. 
Therefore, we can continue (114p by 

$(e,0) < (- (2 + -^—) +l) supVminil^a^il^-^mmilJgV 1 ^!! -27 }- 
Now frrT]) follows from Lemma I3.2( the proposition is proved. 



□ 



The strengthened estimate in the special case of a shearing parameter sequence {sk}k with 
Sk = k for k G 1 will become important when estimating the frame bounds of the concrete 
class of compactly supported shearlet frames we will introduce in Section HI 

The next result reveals the dependence of R(c) on the sampling matrix M c and provides 
a useful estimate for such. For its very technical proof, we refer to Subsection 15.1.21 

Proposition 3.3. Let <S"H(A, ip) be a feasible shearlet system with associated sampling vector 
c = (ci,c 2 ) € (M + ) 2 satisfying c\ > c 2; and let R(c) be defined as in (EJ). Then, for any 7', 
which satisfies 1 < 7' < 7 — 2, we have 



R(c) < Ti£>i(7) + min 



Cl 

c 2 



2 ^ T 2 D 2 ( 7 - 7') + T 3 (£>i( 7 ) + 7J 2 ( T )), 



(15) 



where 



T 1 =p[lC(i) 



logj 



T 2 = P pC(V) 2 lo gi ^ 



+ 



/2ciy 

1 - /iT ' 1 - /i Q -^y \g'/ ' 

/2gc2 

1 - /i 7 ' ' 1 - /i Q " 7 ' ' 1 - /i 7 ' 1 - /i Q " 7 V V q'r 



1 1 1 

zr + " r + 



7-7 



and, /or any h > 0, 

W) = 2(1 + 5^)^(1 + j^), = 6(1 + ^+^0 



/i - 2 
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In particular, for any 7', which satisfies 1 < 7' < 7 — 2, there exist positive constants K\ and 
n 2 (in particular, independent on c) such that 

™^(t) 1+ »&F- (16) 

3.2. General Sufficient Condition. Retaining the notions introduced in the previous sub- 
section, we can now formulate a general sufficient condition for the existence of cone-adapted 
irregular shearlet frames. We defer the lengthy proof to Subsection 15.1.11 

Theorem 3.4. Let ST-L(A,tp) be a feasible shearlet system, let Lj n /, L sup be defined as in 
(jBJ), and let R(c) be defined as in (jUJ). If L sup and L in f are chosen such that 

R{c) < Li n f — Linf and L sup < L sup , 
then S7-L(A, ip) is a frame for L 2 {C) with frame bounds A and B satisfying 

^m&« - R{c)] £ A £ B £ udm\ 11 -' + * (c)| - 

The role of the determinant of the sampling matrix M c , which can also be regarded as the 
inverse of the Beurling density of the translation sampling grid in spatial domain, requires 
a careful examination. The 'quality' of a frame can be measured by the magnitude of the 
quotient of the frame bounds B/A, since the closer it is to 1, the better frame reconstruction 
algorithms perform. In this case, this quotient is bounded by 

B < L sup + R(c) 
A-L inf -R(c)' 

We remark that in order to estimate the quotient B/A, it is sufficient to estimate R(c) and 
Lmf i n (|T7j) . since we can use Proposition 13. II for L sup . 
We make the following observations: 

(a) Notice that as c\ approaches - allowed by (1T61) - the quotient B/A approaches 
Lsup/ L inf , which is optimized by choosing L inf = L inf and L sup = L sup . This is very 
intuitive, since as c\ — > the translation grid becomes denser, the frame property 
more likely, and it becomes merely a question depending on the values Li n f and L sup . 

(b) As we indicated in the previous remark, one can achieve better frame bounds by 
simply choosing a sufficiently small sampling constant c\ (= max{ci, C2}). However, 
this leads to a highly redundant system, which causes a significant computational 
complexity for practical applications. Therefore, we might not want to choose c\ too 
small to avoid high redundancy, but instead seek a threshold, which gives the largest 
possible value for c\ and also c 2 . 

(c) Reasoning about the previous two observations it becomes apparent that the balance 
between the possible values of c\ and c 2 is an optimization problem, which presumably 
is different depending on the application at hand. 
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4. Compactly Supported Shearlet Frames 

We must now brace ourselves for the considerably harder challenge to construct compactly 
supported (separable) shearlet frames with provably 'reasonable' frame bounds. However, 
this task seems, at least to us, much more rewarding as it shows the delicacies of exploiting 
the previously derived general sufficient frame conditions. 

The main idea will be to start with a maximally flat low pass filter and use this to build 
a separable shearlet generator, which is then shown to generate a shearlet frame with 'good' 
frame bounds - in the sense of a small ratio. The modeling of the shearlet mimics the 
classical Example 12.31 with a wavelet-like function in horizontal direction and a bump-like 
function in vertical direction (on C), while ensuring spatial compact support. This is first 
performed for L 2 (C), followed by a discussion of L 2 (M 2 ). The new shearlet frame is finally 
shown to provide (almost) optimally sparse approximations for cartoon-like images. 

4.1. Shearlet Frame for L 2 [C). We start by constructing a compactly supported (sepa- 
rable) shearlet frame iS"H(A, ijj) for L 2 (C). Our construction will ensure that SH(A,ip) is 
a feasible shearlet system, which enables us to exploit our previous results, in particular, 
Proposition I3.1[ Theorem I3.4[ and Proposition 13.31 to derive estimates of its frame bounds. 
For our construction, we focus on regular discrete shearlet systems, i.e., we choose 

A = {(2" j , k2 j/2 , S k2J/ 2A 23 M c m) :j>0,\k\< \2 j/2 ] , m £ Z 2 } with c = (a, c 2 ) £ (M + ) 2 , 

hence the main objective will be to construct a compactly supported shearlet ip, which 
satisfies the decay condition fl6]). 

Our first goal towards this aim is to define a particular univariate scaling function, which 
will become a main building block for the compactly supported shearlet ip we aim to con- 
struct, and prove controllable estimates for its function values. For this, given K,L £ Z + , we 
first let m be the trigonometric polynomial with real coefficients which satisfies m (0) = 1 
and 

\m m 2 = (cosKO) 2 * £ ( K ~l + U ) (sinKO) 2 ". (18) 

n=0 ^ ' 

Notice that mo can be obtained from (fl8|) using spectral factorization (see [5]). The two 
parameters K and L determine the level of flatness near | and 0, respectively. In particular, 
if K = L, then |?n (£i)| 2 is symmetric with respect to the point £i = |, and in this case m 
generates the Daubechies scaling function [5]. The more general trigonometric polynomial 
m we are considering here is sometimes also referred to as the maximally flat low pass filter. 
For an illustration of |mo| 2 we refer to Figure [3j 

The following result states some useful properties, which will become handy in the sequel. 
For the lengthy proof, we refer to Subsection 15.2. 1[ 

Lemma 4.1. Retaining the previously introduced notations, we have the following conditions. 

(i) |m | 2 is even. 

(ii) |mo| 2 is decreasing on (0, \). 

(iii) If K T^_ 2 > \, then \m \ 2 is concave on (0,|). 
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Figure 3. Graph of |m | 2 with K > L (solid), K = L (dot), and K < L (dash). 
The next step is to define the sought scaling function <fi G L 2 (M) by 

oo 
3=0 

Our first result investigating this function concerns a lower bound for the function |0| 2 on 

l 6 ' 6 J • 

Proposition 4.2. Retaining the previously introduced notations, if K ^_ 2 > \, then 

10(6) I 2 > n \m (V) f • e-^CH^fJ)!') . X[ _ 1/6m (^), 6 e R, (20) 

where Jq is chosen sufficiently large. Moreover, the function <fi is compactly supported. 

Proof. First note, that it suffices to consider the interval [0, |], since |</>(£i)| 2 is even (see 
Lemma l4~TT i)). Now assume that £i G [0, |]. Since, by Lemma l47TT ii) + (iii) . |m | 2 is monotone 
and concave on (0, |), we have 

|m (6)| 2 > K(0)| 2 - [6(|m (0)| 2 - |m (|)| 2 )] U- 
Since by construction |mo(0)| = 1, this implies 

imoceoi'^i-ecxieii, 

where 

Ci = l-|mo(i)| 2 >0. 

Let now Jo > be chosen sufficiently large such that 2~ Jo C\ < |. Then, since 1 — x > e~ 2x 
for x G [0, |], we can continue our estimation by 

|m (2- Jo ei)| 2 >e- 12 ^ JoCl )^'. 
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By again using Lemma [4.1( ii). this finally implies 

Jo— 1 00 

|fe)| 2 = [ |m (2-^ 1 )| 2 ni^o(2" J °-^ 



)\ 2 



j=0 j'=0 
Jo— 1 „ 00 

^ nh(¥)| n e^w-'-^ 

j=0 j'=0 
Jo-1 

= n 

j=0 

s nh(¥ 

j=0 



m 



2 e -2- J o+ 2 (6Ci)6 



2 e -2- J o+2Ci 



Finally, by standard arguments (cf. [5]), the function <fi is compactly supported. The 
proposition is proved. □ 

In addition, an upper, 'mathematically controllable' bound for the values of |0| 2 will be 
required for the analysis of the sought compactly supported shearlet ip. We start with some 
preparation for being able to state the precise estimate. For this, we need to introduce a 
trigonometric polynomial ffiQ slightly differing from the previously discussed m®. Letting 
K' e Z be such that < K' < K, we set 

\rh m 2 = (cos(vr6)) 2A " £ (* " 1 + ") (^KO) 2 " 

n=0 V n / 

and again can obtain m from this using spectral factorization. The next estimate for the 
values of |mo| 2 will be crucial for our upper estimate for \<p\ 2 . The very technical proof is 
deferred to Subsection 15.2.21 

Lemma 4.3. Assume that L > 6 and L + 1 < K < 3L — 2 in the definition of mo, i.e., in 
(fl8l) . Let K' be a non-negative integer which satisfies j§rjiz\ — \- Retaining the previously 
introduced notations, we have the following conditions. 

(i) |mo| 2 is increasing on (0, |). 

(ii) We have 

maxK({ 1 )P<E( „ )(^) (sq^) =G- PI) 

71=0 

This now enables us to state the upper estimate for |m | 2 we aimed for. 

Proposition 4.4. Retaining the previously introduced notations, and letting the non-negative 
integers K, K', and L be chosen such that L > 6, L + 1 < K < 3L — 2, and I ^7^_ 1 > \, 
then, for J\ sufficiently large, 



|fei)| 2 <min|l,C 2 -|27rei|- 27 - J[ m (^j 



2 



e 2- J i +1 E„IM«)lh k ) ^ e 
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where 7 = (K — K') — \ \og 2 {C2) and h{n) being the Fourier coefficients of |m p 



Proof. First, observe that |0(£i)| 2 — 1 f° r all ^1 G R, since |m (£i)| 2 < 1. Now fix £1 G 
If |27r£i| < 1, we use the easy conclusion from the proof of Lemma l4.3( i). 

c 2 -J[ 

j=0 



m 



2~3 
27 



; 2- J i +1 E„IK")IH > 1 



to show that 

mini |2tt£i| 



J1-1 

2t n 

J=0 



m 



2~3 
~2^ 



^- J i+ l J2 n \h(n)\\n\ 



l>l^l)| S 



which was claimed. 

Now assume that |27r£i| > 1. First, we observe that 



and 



|m (6)r = cos(7r6) 2(M,) |mo(ei)l 



I0(6)| 2 = II (cos(7T2-^i)) 2(X " X,) |m (2-^i)| 2 . 



Using the classical formula YYjLi cos(2 3 x 



3=0 

00 tn— j _\ sin(x) 



(see 0), it follows that 



sin(27r^i 



2tt6 



2(JC-K') 



(22) 



n 1^0(2^1) 1 2 . 

j=0 

Since |27r£i| > 1, there exists a positive integer J such that 2 J ~ 1 < 2tt|£i| < 2 J . Fixing this 
J, by ( 12 2 p and the definition of C2, we have 



J-2 



|fei)| 2 = |2< l r 2 ^')][|m (2-^i)| 2 |™o(2^i 



j'=0 
J-2 



J=J-1 



< |27reir 2( ^' } J] 2l0g2(C2) II l^o(2- i (2- 7+1 ei))| 2 

j'=0 i=0 

00 

< |27r6r 2(K - K,) (2( J - 1 ) lo ^( C2 )) J] |m (2^'(2- J+1 ei))| 2 . 

3=0 

Since 2 J - X < 2tt|^i| < 2 J implies 2( J - 1 ) lo ^ C2 ) < ^vr^l 1 ^ 2 ^ 2 ), letting 771 = 2- J+1 £ 1; hence 



2 ^< |7?i| < J, yields 



|fe)| 2 < ^r 2 ^-^ 10 ^) SU p I|m (2-^)P 

me[o,iA] j=0 



< C 2 -|27r£i| 



-2( M ')+iog 9 (c a ) sup T| 7 r t0 ( 2 -i ?7l )| ; 

ijl6[0,l/2jr] j=Q 
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By Lemma l4.3( i). this implies that, for all J\ > 0, 



Ji-i 



l</>(6)r <c 2 -|2< 

Since 

|™o(£i)| 2 < 1+ ( 2tt^ 

n 

we can continue by concluding that 



m 



2~i 
2^ 



sup [\rh (2 j Jl r/i) 

«H6[0,l/2jr] j=0 



\h(n)\\n\ |a|<e 



;a«-E„|fc(»»)IM)lfi 



|0(6)r<C 2 -|27r6 



Ji-i 

21 n 

3=0 



m 



27 



• e 



2- /l+1 E„IM") 



The proposition is proved. 



□ 



We can conclude from this result that |0(£i)| = 7 ) for < K' < 4K 5 L+1 (which is 

equivalent to R ,^_ x > \), where 7 = (K — K') — | log 2 (C2). Let us now take a closer look 
at the decay rate 7. This depends on K' for fixed K and L, where K' can be chosen to be 
an arbitrary integer satisfying < K' < K. Now regarding the parameters K' and K as 
control parameters for the low pass filter mo given by 



L-l 



|m (£i)| 2 = cos «i 



XK-K') 



COs(7T^i) 2 ^ ^ 



n=0 



K 



-1 + n 



(sin«i)) 



2/i 



we see that, for any choice of K' satisfying < K' < K, the same scaling function <fi is 
generated. Therefore, for analyzing the decay of <f>, it is sufficient to show that there exists 
some K' such that 7 is above a certain threshold; later we aim for 7 > 3. The following 
lemma makes this consideration explicit by choosing K' — 0. We however would like to 
emphasize that this choice was made only for proving feasibility, i.e., ([6]), of the shearlet 
which we will introduce in Theorem 14. 7j it will not be our optimal choice when estimating 
L 

to Subsection 



;p and R(c) of the generated shearlet frame. The technical proof of this lemma is deferred 



Lemma 4.5. Let trq be the low pass filter defined in fTLSl) with K > ^ and L > 2, and let 
C2 and fho be as in Lemma \4-3\ with K' = 0. Th 



en 



max |m (£i)| 2 < 2 2R - L ' 2 - 1 
£ie[0,i] 



and the constant 7 from Proposition \4-4\ satisfies. 

7 



K- 1 -\o g2 (C 2 )> 1 -(^ 



Finally, we have reached the stage, where we can introduce a compactly supported shearlet 
ip G L 2 (R 2 ), which generates a (separable) feasible shearlet frame iS"H(A,^) with 'good' 
mathematical controllable frame bounds. For illustrative purposes, some elements of the 
shearlet frame iS"H(A,?/;) are displayed in Figure HI 
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Proposition 4.6. Let K, L G Z + be such that L > 10 and ^ < K < 3L — 2, let m be the 

associated low pass filter as defined in (I18p . and let (j) be the associated scaling function as 
defined in f|T9|) . Further, we define the bandpass filter mi by 



|™i(6 



|m (ei + l/2)| 2 , 



Then the shearlet defined by 

m = m L (46)<fe)!(26): C = (6,6) e K 2 
satisfies the feasibility condition ([6]); m particular, for any < K' < 4K ~ L+1 ; we have 



where 



|^(OI < min{l, | g ei| a }min{l, |g'6|" 7 } min{l, |r£ 2 p}, 



a = K — K' and 7 = (iT - if') - - log 2 (C 2 ) 



(23) 
(24) 



wift C2 fremg defined in (|2ip . and 



J1-1 



g = 4vrC 2 2(K ^ , ,9 / = 27r(C 2 n 

3=0 



m 



2-J 
2^T 



2n 



; 2- J i +1 EnW«)IN ; r = 2g '. (25) 



Further, we have 



'Jo-l 



m)\ 2 >\mo{\)\ 2 - [ n 

.3=0 



m 



2^ 

IT 



-2-^0+2(l-| mo (l)|2) 



e 

^erefi = {£=(6,6)eM 2 : 6 G [-L ±] U 6 G £]} 



Xn(0 > 0, £ G 



(26) 





1 


\\ 


1 


II 


1 


// 


1 


✓ 



Figure 4. Display of compactly supported shearlets ip2,k,o, k = —2, —1, 0, 1, 2, 
in spatial domain. 

Now for K,L G Z + satisfying L > 10 and ^ < K < 3L — 2, Proposition 14. 61 together with 
([TT]) implies that L sup < 00. Moreover, by Proposition 13.31 -R( c ) can be made arbitrarily 
small by choosing a sampling matrix M c with sufficiently small determinant. Finally, let the 
set Q be defined as in Proposition 14.61 From ( 1261) . we see that |^(£)| 2 > -^m/ ' Xn(£)> where 

/Jo-l 

^m/ := |m (|)| 2 JJ 

\3=0 



m 



_2-^0+2(l-| mo (I)|2) 



(27) 
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Since 

oo 

[J |J A 21 SlQ = C, 

3=0 |fe|<[2J/ 2 ] 

it then follows that 

mO) = J2\^(SlA 2 ^ ! :)\ 2 >L mf -J2xn(SlA 2 -^)>L mf on C. 

j,k j,k 

Thus, by Theorem 13.41 both the lower and upper frame bound exist with explicit estimates 
for the shear let frame iS"H(A, ip) for a sampling matrix M c with sufficiently small determinant. 
This directly implies the following main result. 

Theorem 4.7. Let 

A = {(2" j , k2 j/2 , S k2j /2A 2i M c m) :j>0,\k\< [2 j/2 ] , m £ Z 2 } 

be the regular sampling set of E> cone . Further, let K,L £ Z + be such that L > 10 and 
^ < K <3L - 2, and define a shearlet ip £ L 2 (R 2 ) by 

^(0 = 7^(4^)0(6)0(26), £ = (6,6) e K 2 , 
where m is the low pass filter satisfying 

|m (6)| 2 = (cos(7r6)) 2 ^ ( K ~l + U ) ( sin «0) 2n > 6 e R, 

n=0 ^ ' 

mi zs i/ie associated bandpass filter defined by 

K(6)l 2 = |m (6 + i/2)| 2 , 6eK, 

and <p is the scaling function given by 

oo 

fe)=n>°( 2 ~^)' ^ gm - 

i=o 

Then there exists a sampling constant c,\ > siic/i that the shearlet system Sl-L{K,ip) forms 
a frame for L 2 (C) for any sampling matrix M c with c = (ci,c 2 ) £ (M + ) 2 and c 2 < c\ < c\. 
Furthermore, the corresponding frame bounds A and B satisfy 

, - TTT^ iLinf - R(c)\ < A < B < , - )„.. [L aup + R(c)}, 
|det(M c )| L ; K JS ~ ~ ~ |det(M c )| L p WJ 

where R(c) < L inf < L inf and L sup < L sup . 

To illustrate how this result can be applied, we discuss the derived estimates for the frame 
bounds using the following particular choices for the free parameters K, K' , L,c\, and c 2 . 
However, we do not claim that this is the optimal choice. In fact, optimizing the estimate 
for the ratio of the frame bounds A/B is a highly delicate task in this situation. 
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K 


L 


Cl 


c 2 


K' 


B/A 


K 


L 


Cl 


c 2 


K' 


B/A 


39 


18 


1.00 


0.40 


(27,17) 


37.1204 


39 


19 


0.9 


0.40 


(27,18) 


44.5359 


39 


18 


1.00 


0.30 


(27,15) 


32.0208 


39 


19 


0.9 


0.30 


(27,16) 


28.4307 


39 


18 


1.00 


0.25 


(27,15) 


31.9105 


39 


19 


0.9 


0.25 


(27,15) 


28.0983 


39 


18 


1.00 


0.15 


(27,15) 


31.9019 


39 


19 


0.9 


0.20 


(27,15) 


28.0699 


39 


18 


1.00 


0.10 


(27,15) 


31.9019 


39 


19 


0.9 


0.15 


(27,15) 


28.0683 



(a) (b) 



Table 1 . Some estimates for the ratio of the frame bounds for the compactly 
supported shearlet frame constructed in Theorem 14 .71 for various choices of the 
parameters K, K', L and the sampling constants Ci, c 2 . K' is given in the form 
(■,■), where the first component is an estimate for L sup chosen by (fTTl) and the 
second component is an estimate for R(c) chosen by (TL5I) . 



Example 4.8. For our discussion, we decided upon the following choices for the various 
parameters necessary for exploiting Theorem \4- 7| to derive an estimate for the ratio of the 
frame bounds: 

We may choose Li n j as in (I27|) for the expression for the lower frame bound. The value 
L sup can be estimated using (JTTJ) with \i = \,p = 1. Further, for a given sampling matrix M c 
with c = (ci,c 2 ) G (M + ) 2 and c\ > c 2 > 0, the value R(c) can be estimated using (Tl5l) with 
\i — ~, p = 1 from Proposition \3.3l For both those estimates, given positive integers K,K', 
and L such that for K' < 4K ~^ L+l we choose \i=\, the values q,q',r are chosen as in (|25|) . 
and the values a, 7 as in (124j) . 

Table \4T8] shows some estimates for the ratio of the frame bounds ofS7i(A,ip) using Theo- 
rem \4- r l\ for specific choices of the parameters K,K', and L with various sampling constants 
ci,c 2 . It should be pointed out that numerical experiments show much better results, but the 
presented ones are the estimates we are able to prove theoretically. Such estimates do not 
even exist for compactly supported wavelet frames! 

Notice the flexibility in choosing K,K', and L. A full blown optimization analysis might 
lead to an even better estimate for the ratio of the frame bounds of SH(A,tp), but this is 
beyond the scope of this paper. 

4.2. Shearlet Frame for L 2 (R 2 ). We now aim to generate a frame for the whole space 
L 2 (R 2 ) using the compactly supported (separable) shearlet frame introduced in the previous 
subsection. One classical way for achieving this is by orthogonally projecting each shearlet 
element in Fourier domain onto the cone C, and proceeding likewise for the vertical cone. 
However, this procedure is quite counterproductive in the sense that it destroys most of 
the advantageous properties, e.g., compact support and regularity, of the elements of such 
a shearlet frame, in particular, if the Fourier transform of the to be analyzed signal is not 
entirely supported in either the horizontal or the vertical cone. 

This problem can though be quite easily resolved by taking the union of two cone-adapted 
shearlet systems - one for the horizontal cone C = C\ U C3 and the other for the vertical cone 
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C = C2 U C4. In fact, a shearlet frame S7i(c;(f),x/),x/)) (cf. Definition 12.21) for L 2 (R 2 ) can be 
constructed by the following 

Theorem 4.9. Let ip G L 2 (R 2 ) be the shearlet with associated scaling function G L 2 (IR) 
both introduced in Theorem \4- 7[ and sei 0(a;i,X2) = 0(xi)^(a;2) and ?/;(£!, x 2 ) = ^(a^,^)- 
Then the corresponding shearlet system SH(c; <p,ip,ip) forms a frame for L 2 (M 2 ) for any 
sampling matrices M c and M c with c = (ci, C2) G (M + ) 2 and c% < C\ < C\. 



Proof. In the sequel, the constants 0,7 and q,q',r will be those defined in ( 1241) and ( 125 
respectively. 

First note that the function $ defined in ([7]) now becomes $ : R 2 x I 2 4 1 defined by 



$(£, w) = |fe) + w)| + $!(£, w) + $ 2 (£, w), (2* 



where 



and 



i>o |fe|<r23/2] 



$ 2 (6^) = E E |^(^2-.0 ^(5 fc i 2 -^ + 0. 

Also, -R(c) defined in ([2]), is now given by the term 

(r (^ 1 m)r (-cr 1 m))' + (r 1 (M c - 1 m)r 1 (-M c - 1 m))^ + (r 2 (M < r 1 m)r 2 (-M ( r 1 m))i 



r (w) = ess sup |fe)| Ife + w ) I an d T^o;) = ess sup w) for i = 1,2. 



where 



In fact, Theorem 13.41 can be easily extended to this case so that we only need to estimate 
Linf by some positive L 

inf 1 L SU p by some finite L sup , and R(c) to derive finiteness of the 
estimate for the ratio of the frame bounds in ffTTl) . Using the function from (I28p . 

Lj„j and L sup are defined as in (JBJ) and ([9]) - note that in this case, C is replaced by M 2 in 
the definition. 

We start by estimating the first term |fe)||fe + u))\ from WLOG, we may assume 
that ||77i||oc — m i 7^ for m G Z 2 \{0}. Then, by Proposition 14. 4^ we have 

Ifejllfe+cr 1 ™)! < 10(6)110(6 + ^)1 

< min{l, |g'6|~ 7 } min{l, \q'^\ + q c^mi^ 1 } 

< feVlHI^ (29) 



q 

provided that c^q'mi > 4. Now choose c\ > sufficiently small so that c~{ q 1 > 4. 

For the second and third term in 02 8|) . we can apply the estimates ffTTl) and ffTB]) from the 
proof of Proposition 13.31 In particular for the third term, we only need to switch the roles 
of the variables £1 and £2 to reach the same estimates as (flTT) and (TIE]) . 
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Set now 7" = 7 — 7' for an arbitrarily fixed 7' satisfying 1 < 7' < 7 — 2. Then, by (J29 
and ( ITEj) . we obtain 



R(c)<C(=^) +C 



2gc 2 



q'r 



(30) 



Also, using ( flTT) . 



Lsup ^ 1 



27 - 1 



log 2 



)l-2a 



+ 



1 - 2- 2 ^ 



(31) 



since < 1. This gives the upper bound of L sup and we may choose the upper bound 



(32) 



L sup to be the right hand side of (j3T|) . 

Finally, let us estimate the lower bound of By (1201) and ( 1261) . we have 



|0(£)| 2 >Cxn o (£), l^(£)| 2 >CW£), and |^(0| 2 > C X n 2 (£), 



where 



^0 
fix 

n 2 



[2 

d2 



iieiu< j>, 

12 < l&l < 6' N < n)' 



k<l&l<i,lfcl<i*}. 



12 



Setting 



we observe that Q 



Q = O U |J Ajy-^A U |J i 2J 5jfi 2 

^ 2 and, by ([32]), 

$(£, 0) > L inf ■ xn{€) fo r some L inf > 0. 



(33) 



This implies $(£, 0) > > on IR 2 , which yields the lower bound of Li n f. 

Concluding, the estimates (1501) . (l3Tj) . and (1331) provide all required constants L in f, L sup 
and the upper bound of R(c) in (fT7j) . Especially, we can choose a sampling matrix M c with 
sufficiently small determinant such that Lj n / — -R(c) > 0. This proves the existence of the 
frame bounds. 

Finally, it is obvious that all functions in ST-i(A, A,ijj,'ijj,9) are compactly supported in 
spatial domain, which finished the proof. □ 

In the proof of Theorem I4.9j we see that our upper bounds of L sup and R(c) are about 
twice as large as the upper bounds in the cone case and the lower bound L in f is about the 
same compared to the cone case. Also c can be chosen so that R(c) is sufficiently small. This 
indicates that in this case, our estimate for the ratio of the frame bounds is about twice as 
large as the estimate for the ratio B/A'm the cone case for sufficiently small determinant of 
the sampling matrix M c . 
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4.3. Sparse Approximation using Compactly Supported Shearlets. One essential 
performance criterion for a frame composed of anisotropic elements - besides the ratio of 
the frame bounds - is the approximation rate of curvilinear objects. This viewpoint arose 
due to the fact that edges are the most prominent features in images, hence representation 
systems should in particular provide sparse representations for those. 

To stand on solid ground, we first briefly recall the mathematical model of a cartoon-like 
image introduced in pQ. For v > 0, the set STAR 2 (v) is defined to be the set of all B C [0, l] 2 
such that B is a translate of a set {x G M 2 : \x\ < p(6),x = (\x\,9) in polar coordinates} 
which satisfies sup \p (9)\ < v, p < p < 1- Then, £ 2 {y) denotes the set of functions 
/ G L 2 (R 2 ) of the form 

/ = fo + fxXB, 

where fo, fi G Cg([0, l] 2 ) and B G STAR 2 {y). In [7], it was proven that the optimal ap- 
proximation rate for such cartoon-like image models which can be achieved under some 
restrictions on the representation system as well as on the selection procedure of the approx- 
imating coefficients is 

\\f-f N \\t<C-N- 2 asiV-^oo, 

where /jv is the iV-term approximation generated by the N largest coefficients in magnitude. 

In [TB], two of the authors proved that a large class of shearlet frames - including, in 
particular, a significant set of compactly supported shearlet frames - provide optimally 
sparse approximations of such cartoon-like images. For the convenience of the reader, we 
state the result below. Notice that in [18], the result was proven for a isotropic sampling 
matrix, but the extension to our anisotropic sampling is immediate. 

Theorem 4.10 ([IB])- Let c > 0, and let (f>,ip,ip G L 2 (M 2 ) be compactly supported. Suppose 
that, in addition, for all £ = (^1,^2) £ K 2 , the shearlet if) satisfies 

(i) h?(£)l < C ■ min{l, • min{l, ■ min{l, |£ 2 |" 7 } and 



0(o\<\hm-{i+ [ § l 



where a > 5, 7 > 4, /i e L 1 (M), and C is a constant, and suppose that the shearlet if; satisfies 
(i) and (ii) with the roles of £1 and £2 reversed. Further, suppose that SH(c; forms 
a frame for L 2 (R 2 ). 

Then, for any v > 0, the shearlet frame Sl-L{c] </>, if), i/j) provides (almost) optimally sparse 
approximations of functions f G S 2 (u), i.e., there exists some C > such that 

11/ - f N \\l < C ■ (log iV) 3 • N~ 2 as N 00, 

where f^ is the nonlinear N-term approximation obtained by choosing the N largest shearlet 
coefficients of f. 

Using this result, we can prove that in fact the compactly supported shearlet frames 
we introduced in Theorem 14.71 even provide (almost) optimally sparse approximations of 
cartoon-like images, i.e., functions in S 2 . Although the optimal rate is not achieved, the 
log-factor is typically considered negligible compared to the A r_2 -factor, wherefore the term 
'almost optimal' has been adopted into the language. 
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Theorem 4.11. Let ip G L 2 (R 2 ) be the shearlet with associated scaling function G L 2 (M) 
both introduced in Theorem \4-7\ with L > 18, and set <p(xi,X2) = <f)(xi)(f)(x2) as well as 
ifj (xi,X2) = ip(x2,Xi). Then the compactly supported shearlet frame iS"H(c, c; <j),ip,ip) provides 
(almost) optimally sparse approximations of functions f G & 2 {y), i-c, there exists some 
C > such that 

11/ -f N \\t<C- (log iV) 3 • N~ 2 asN^oo, 
where J'n is the nonlinear N-term approximation obtained by choosing the N largest shearlet 
coefficients of f . 

Proof. The fact that ST-L(c; <f),ip,ip) is a compactly supported frame follows from Theorem 
14.91 Hence it remains to prove that the shearlet if; satisfies conditions (i) and (ii) from 
Theorem 14.101 and the shearlet ip likewise. 

By Lemma 14.51 and Proposition 14. 6[ there exists a constant C such that 

|^(0| < C • min{l, U K ] min{l, \^\- L/ '~ 1/2 } min{l, U^ 1/2 }, 

where L/4 + 1/2 > 4 and K > 5. This already proves condition (i) of Theorem 14.101 

To show condition (ii), we first observe that there exists a function ipi such that we can 
write = ^1(^1)0(2^2) so that 

|^i (6)| <C-min{l,|ei|^}min{l,|eir L/4 " 1/2 } and |0(£ 2 )| < C • min{l, U^^ 2 }. 

Hence it is sufficient to prove that 

|(#)'(6)| < C ■ min{l, |6|- 7 '} with i > 4. (34) 

For this, we first choose a positive integer p' > such that 4 < p' < -| — ^ < K, which we 
are allowed to do since L > 18. Then, by the definition of 0, 




where |m (6)| 2 = En=o ( K ~n +n ) (sin« 2 )) 2n . Now define S and by 

/sim>6)y' ~ / sm(27r5) \ K ~ p ' _ 

S ^)={-2^r) and ^ 2)= {-2^2~) n-o(2^ 2 ), 

which gives 

(0)'(6) = 5 / (6)0(6) + 5(6)(0)'(6)- 

Since max{|S"(^2)|, l-S^)!} < C\(,2\~ p , it remains to show that both |0| and |(0)| are 
bounded. Notice that we may assume that m (0) = 1 and is compactly supported (see 

[5]). Therefore, boundedness of implies that and (0)' are bounded, which enables us to 

restrict our task further and it remains to check that G L X (1R). This will now be proved 
by showing that, for all £2, 

10(6)1 < C-mm{l,\^ v } with 77 >1, (35) 

which obviously implies G L X (R). 
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First, it is easy to show that \<ft\ is bounded on [—1,1]. Therefore, if |^| < 1, then fl35|) 
holds. Now assume that > 1- Then there exists a positive integer J > such that 
2 J_1 < < 2 J , and, using the fact that Lemma I4T51 implies max^ 2 \fho($,2)\ 2 < 2 2K ~ L / 2 ~ 1 , 
we obtain 



j-i 



10(6)1 < ci^i-^n^ximo^ini^ocs-^ 

?2 



j=0 



j=0 

00 



< C| G |-(^-p'-W + i/2) sup J|^ o(2 -^ 2 



Ce[-i,i] 



3=0 



-(L/4+1/2-p') 



< c|6| 

Hence, L/4+1/2 — p' > 1 and this proves (135 p . and hence (I34p . The theorem is proved. □ 

5. PROOFS 
5.1. Proofs of Results from Section [31 

5.1.1. Proof of Theorem \3.4\ Let / e L 2 (C). Then, by definition of ipj t k,m, 



E E E | (/' fem>i 2 (c) 

jgz fceXj mez 2 



EE E 

j>0 fcGXj meZ 2 



3/2 



(36) 



We now first decompose the sum over m G Z 2 . For this, set ft = [-5, 5] . Then, by 
appropriate changes of variables, 



E 



3/2 



E 



E 

m€Z 2 

E 

mez 2 



3/2 



det(M c 
-3/2 



/ ( a-^s^m;^ ) X c ( A- L s^'M^ ) v (M-^y m ^ m} d£ 



-1 0— T A/r— li 

»3 Sfc 



-lCU27ri($,m> 



det(M c 



3/2 



det(M c 



E / (^fX^ + £)) xc (a^S~ t M-\Z + , 
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Next, applying Plancherel's theorem, 



E 

mGZ 2 



3/2 



c 

3/2 



|det(M c )|y n 



E / ( A a-S7 k T M^ + £)) xc (a£S? M- 1 ^ + £)) ^ (M-i(e + £)) 

fez 2 



Resolving the absolute values yields, 



E 



3/2 



-3/2 

iJRl/ ^(W^e + fl)»G 



•/ (A-^-fM-He + m) J ^ + m)) 

-3/2 



|det(M c )| 



E / (A£S 7 TM-i(t + m-£))$ (M- 1 ^ + m - £)) 



-3/2 



|det(M c )| 7 R2 



E / (a£s£M?£) xc (A-'S^ M-^) / (A-^-^M-He + m)) 

meZ 2 



Combining with (136|) . 



E E E |(/'^^)^(c) 

jeZ fceiC, mGZ 2 



|det(M c )| 



£ jez fceKj mez 2 



•i> {Sf k A aj ^ {SjA a ^ + M~ x m) di 



T x + T- 



2- 



(37) 



where 



^^^EE/NV^) 
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and 



det(M, 



m 



j&L k&Kj ° c mGZ 2 \{0} 



By using Cauchy-Schwartz inequality, 

1 " s ' 12 E [rlM^r^M-V)] 172 

meZ 2 \{0} 



IT, I < 



det(M c ; 

By (1571) and (1551) . we finally obtain 



(CO 



jSZ fceKj m.gZ 2 



< 



L S up ~~ -R(c) 

|det(M c )| 



The introduction of the lower bound for Lj n / in the estimate 



immediate. The theorem is proved. 



det(M c )| 



(38) 



[L inf - R(c)] is 



5.1.2. Proof of Proposition [PI We start by estimating T(2ui, 2w 2 ), and will use this later 
to derive the claimed upper estimate for R(c). For each (wi,u; 2 ) € IR 2 \{0}, we first split the 
sum over j as 

T(2lu u 2u 2 ) 

= ess sup E E r s k a j£i + a i /2 ^) i> (a>j£i + 2wi, s^-fi + a 1 / 2 ^ + 2^ 2 



{j ; K?i|<IMIoo} {j:K?i|>IMIoo} / fcsz 

l/2 ( 



V> + 2wi, s fc a^i + a/ £ 2 + 2w 2 
ess sup (Ji + 7 2 ), 



(39) 



where 

Ji = E ^ ( a i£i> s fc a i£i + a i /2 ^) ^ + 2wi, s fc a^i + a] /2 ^ 2 + 2w 2 
and 



■l + a) /2 6 + 2cj 2 



O':|a^i|>l|w||oo} ^62 

The next step consists in estimating I\ and J 2 . Before we delve in the estimations, we 
first introduce some useful inequalities which will be used later. Recall that a > 7 > 3, and 
q, q', r are positive constants satisfying q > q' > 0, q > r > 0. Further, let 7" = 7 — 7' for 
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an arbitrarily fixed 7' satisfying 1 < 7' < 7 — 2. Then we have the following inequalities for 
x,y, z G R. The proofs are all elementary, and hence we will skip them. 

min{l, |gx| Q }min{l, |g'x|~ 7 } min{l, |ry|~ 7 } 

< min{l, |gz| Q_7 } min{l, |g'x|~ 7 } min {l, |(gx)~ 1 ry|" 7 } , (40) 



min{l, \x\ 7 } min < 1 





1- 


-z 














x - 







< 2 7 "|y|- 7 "min{l,|x|- 7 '}max{l,|l + zr"}, (41) 
min{l, |gxr~ 7 } min{l, |g'xr 7 }|x| 7 " < (g')~ 7 " , (42) 



and 



min{l, |gx| Q " 7 } min{l, |g'x|- 7 }|x| 7 " < (g')~ 7 " min{l, \qx\ a ~ 1+1 "} min{l, |g'x|~ 7 '}. (43) 
We start with I x . By applying §6§ and ( l4"Uj) . 

/x < ^ min{|ga^ir 7 , 1} min{|g'a^i|- 7 , 1} ■ min{|g(a^ 1 + 2^)| a " 7 , 1} 

{j:K?l|<IMI°o} 



■ min{|g'(aj£i + 2cdi)| 7 , 1} y^min 



6 



1/2 e 



■ mm 



2a; 2 



6 



1/2 1 



1 



2wi 



1 



(44) 



We now distinguish two cases, namely ||w||oo = |^i| > an d IMloo = l w 2| > l a jCi|- 

Notice that these two cases indeed encompass all possible relations between u and £1. 

Case 1: We assume that \\loWoo = > |aj£x|, hence |a 3 -£i + 2tui| > |a;i|. This implies 
\q'{a£i + 2^i)p < ||g'w||~ 7 . Thus, continuing (ill), 

Jx < ]T min{|ga J £ i r- 7 ,l}min{|gVir 7 ) 1 }k'(^i + 2 ^i)r 7 



{j : |aj£l|<IMI°o} 

■E 



mm 



£2 



1/2/. 



< H^ll" 7 minllga^l^alminllgVil" 7 ^} 



{j:|ajfi|<||w||oo} 

q \ - r 



■ - > - mm 



r 14 — ' (/ 

fcez * 



r 

- I s k 

q 



1/2 1 



n. 



By (jlj) and Lemma [3.21 finally conclude 

(^(7)) 



h < 



\\q'w\\L 



J^mindga^il" 7 , 1} min{\q'a^x\ 7 , 1} 



lly ^lloo 



log 



1 1 

+ z + 



l-ff 1 - n a - 



(45) 
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Case 2: We now assume that H^Hoo = (cc^l > Then (jHJ) applied to 0HJ yields 

h < 2 7 " min{|ga^ir-^l}min{|gVil" 7 ; l}min{|g(a : ,ei + 2u; 1 )| Q - 7 ,l} 

l?:|aj£i|<IM|oo} 



min{|g'(a^i + 2oj 1 )\ 7 , 1} y^min 



fcez 



- s k 



6 



1/2 1 

a/ £i 



2u 2 


-V , 






2wi 






max < 


[•■ 


H 















Applying now (j4j), we obtain 

(W)) 



Ji < 2 7 



f ^min-^a^iP 7 , 1} min{|g'a^i| 7 , 1} min{|g(a i ^ 1 + 2wi)| a 7 , 1} 



■ min{|g'(aj£i + 2ux)\ 7 , 1} |a,-£i| 7 max jl, 

We next split Case 2 further into the following two subcases. 
Subcase 2a: If 1< 11 + ^H, then 



2wi 



(46) 



max < 1, 



(■■ 








«i6 





— I a j£i + 2wi | 



Hence, by Lemma [3.21 and exploring inequality (|42p . we can conclude from ( 146]) that 



(aery)) 



— W||oo 



log 



l//x 



1 1 

+ 



1 - 1 - fi a ~ 



(47) 



Stz&case 2b: If 1 > |1 + then, for all j > 0, 

min{|g(%£i + 2wi)| a ~ 7 , 1} min{|g'(a^i + 2wi)|~ 7 , l}max ^ 1, 



!■ 




7 "1 




Ufa 





< 1. 



Hence, by exploring inequality (S3]), we can conclude from (T4"6"]) that 



r < (?g(V)) 



9 J", ,117 



I? 



7 + 



1 - /iT' 1 - / ' 



(48) 



We next estimate /2- First, notice that the inequality (1441) still holds for /2 with modified 
index set for j - we have {j : |aj£i| > IMloo} instead of {j : |aj-£i| < ||w||oo} in this case. 
Therefore, by (EQ), 

h < 



^2 minllga^xl 7 ,l}min{|g / a j ^ 1 | T ,l}^min< 

{j:|«j^i|>II<^IIoo} fcez L 



Sk + 



1/2 j. 

a/ £i 
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^ il C ^) E min{|ga^i| Q - 7 ,l}min{|g'a^ir 7 ,l} 

{j:K£l|>||w||oo} 

< (^w) E i^Vir 



{j' : l a i?i|>ll' il lloo} 



By Lemma 13.21 this estimation can be finalized to 



h < Pipy 



g'cu||2o VI - ffi 



(49) 



We are now ready to prove the claimed estimate for R(c) using ( 1391) . For this, we define 
the constants T 1; T 2 , and T 3 as in Proposition 13.31 Now define 

Q = {m G Z 2 : \m\\ > |^ 2 |} and Q = {m G Z 2 : |mi| > c^" |m 2 |}. 

If c^ l \mi\ > c^ 1 \m 2 \, then, by (jUJ) and (1391 . 

r(±M c _1 m) < TilHI^ 7 + T 3 ||m||~ 7 for all m G Q. 

If c^lmxl < c 2 - 1 |m 2 | with m ^ 0, then, by (l47lh (j48|), and (l49l. 

T(±M^ 1 m) < T 2 \\m\\^<" + T 3 ||m||- 7 for all m G Q c \{0}. 

Therefore, we obtain 



R(c) 



(r(M- 1 m)r(-M- 1 m)) : 

m€Z 2 \{0} 



< I J^Ti||m|| 00 7 + T 3 ||m| 



-7 
oo 



E T 2lHloo 7 " +T 3 ||m| 

^mGQ c \{0} 



Notice that, since Q C Q, 



Also, we have 



E 



|m|| 7 < ||m| 



||m|| 7 < min 

meQ c \{0} 

Therefore, (1501 can be continued by 



£i 

^2 



]T ||m| 

meQ c \{0} 



R(c)<T 3 



\m\ 



1 7Ti 1 1 + min 



meZ 2 \{0} 



£l 

c 2 



2 71 2 



mGS c \{0} 



(50) 
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To provide an explicit estimate for the upper bound of R(c), we compute XlmeS ll m lloo 7 an d 
E m gQc\ IHL 7 as follows: 



and 



oo oo oo 

\mi\ 

m£Q mi=l m2=lmi=m2+l 



1 \ 4 / 1 
< 2 1 + -)+ -II 



7 — 1/ 7 — IV 7 — 2 



£ Nil* = Eimi« +4 £ 

meS c \{0} meQ mi=l 



oo 

'mil 



1 \ 4 / 1 

6 1 + r)+ 1 



7 — 1/ 7 — 1 \ 7~2 
This completes the proof. 



5.2. Proofs of Results from Section 14.11 



5.2.1. Proof of Lemma 1^.1, First, it is obvious that the function mo is even. Next, letting 



y = sin 2 (7r£i), the values of the trigonometric polynomial |mo(£i)| 2 can be expressed in the 
form 



m 



(6)| a = P(y), where P(y) = (1 - y)*^ ( K 1 + ") 

n=0 ^ / 



We compute 



P'iv) 



(i-yf- 1 

(i-yf- 1 
(i - y) K ~ l 



K-l + n 



n 



/ n + (i-y)E 



K - 1 + n 



n 



,n-l 



L-l , Tr \ L-l 

n=0 v n=l 

T ^ fK-l + n\ n ^ (K + n-l\ n 

A £ n )y -E n ( „ )y B +E(«+i)( . , i. 

n=o v / n= i v / n= o 



71+1 



—K 



K + L-2\ L _ x .(K + L-2 X r_ x 



L- 1 



V 



(L-l) 



L - 1 
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Hence, P'(y) < for y £ (0, 1), and this immediately implies that |m (£i)| 2 is decreasing on 
(0, |). The second derivative can be derived by 

P"{y) 

^ ■^l_ y) K-i + {K + L _ l){K _ l] 



-(L-1)(K + L- 1) 



y 



K + L-2 
L—l 

(K -1)(K + L-1 



L - 1 

y L -\i-y) K - 2 

K + L-2\ (K + L-2 



K-l 



L — l 



y L - 2 {i-y) 



K-2 



L—l 



K + L-2 



Thus P"(y) < for < y < j^j^. Hence 



d 2 , 1 

rrzoC^i) | 2 < for < £i < — arcsin 

71 



Since 



L-l 
K+L-2 



> \ implies 



L-l 



K + L-2 



arcsin 



7T 



L-l 



K + L-2 



> 



6' 



concavity of |mo| 2 on (0, h) is proven. 



5.2.2. Proof of Lemma \4-3\ Letting y = sin 2 (7r£i) for £i £ [0, 1], the values of the trigono- 
metric polynomial |m (£i)| 2 can be expressed in the form 



L-l 



\mo(Hi)\ 2 = P(y), where P(y) = (1 - y) K ' ^ 



n=0 



if - 1 + n 



n 



V 



If if = 0, P(y) is increasing on (0, 1), which proves one part of claim (i). 

To prove the remaining part of claim (i), let us assume that K' > 0. For y £ (0, 1), direct 
computation yields 

(P)'(y) 



L-2 



[l-yf- 1 (K-K')^: 



n=0 



K-1 + n 



n 



y n -{K' + L-V 



K -L-2 
L - 1 



.L-l 



y 



L-2 



-\i-yr~ 1 (k-k')y: 



n=0 



K — 1 + n 



n 



,n-L+l 



(K' + L — l) 



K -L-2 
L- 1 



Hence, (P)'(y) is a product of a decreasing function and of y L ~ l (l — y)^'" 1 > on (0, 1). 
This implies that P(y) is increasing on (0, z ) as long as (P)'(z ) > 0. Obviously, the same 
is true for \m (^i)\ 2 . Thus, since ~ = sin 2 (|), it suffices to show that (-P)'(j) > 0. For this, 
we compute 
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Since A > 7 and K \ L ~[ 2 < 4, the term in the bracket [ • ] can be further estimated by 
22L f K-K> ^ Q g (K - 1 + n\ 2 „ 2n /A' + L - 2 



K' + L- 1) \AJ ^ \ n / L-l 

n=0 



2L 



> 2 



> 2 2i 



A' — K 



K' + L 



'\ (l\ K ( K + l ) , A-(A+l)(A + 2) 

■l/Wi 4 4 2 -2 4 3 -6 



A + L - 2 
L - 1 



A -A' \ (K + L-2 



K' + L-lJ V L-l 



K ~ K ' \ - + > W ) _ 2— > 0, 



,A' + L-1/ (L-iy. - VA' + L-l, 

where the last inequality follows from ^tF^zi — I" This completes the proof of claim (i). 



To prove claim (ii), observe that, for < y < 1 



n=0 
L-l 



n=0 

Thus claim (ii) is proven. 



A - 1 + n\ ( K' \ k' / n 



n )\K' + nJ \K' + n 



5.2.3. Proof of Lemma \4 . 5\ Without loss of generality we assume that L is even, since the 
'odd case' can be proven similarly. Since A' = and by obvious rules for binomial coefficients, 

n=0 v 7 n=0 v 7 n=L v 



Utilizing the estimate J2n=o ( K ~l +U ) < 22K ~ 2 from 0' we have 



vR"-l (K-l+n 

n 



ft<2»--E(* r " 1+ 1 <«) 



n=L 



Now let a n = ( A ' ^ +n ) for n = 0, . . . , L — 1, and note that 

a-n+i _ A + n > K + L-l > 
a n n + 1 ~ L 

which implies 

dL + --- + a>K-i > a L + a L+ i + ■ ■ ■ + asi, 

2 x 

> 2 L " 2 (2ai) + 2 L " 3 (2a 3 ) + • • • + 2 L/2 -\2a L . l ) 

> 2 L ' 2 - 1 (a + --- + a L - 1 ). 



Thus, using f[5"Tj) . 



C 2 < 2 2A '" 2 - 2 L / 2 - 1 C 2 
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and hence 

q < 2 2K ~ 2 - L / 2+1 — 2 2K ~ L / 2 ~ 1 

Thus, we have 



max|m (6)r<2 2i "" i/2 " 1 
and 

2 7 = 2K-log 2 (C 2 )>~ + 1, 

which is what was claimed. 

5.2.4. Proof of Proposition \4~J^ By (i) and (ii) in Lemma [4.11 we obtain 

|mi(4£i)| 2 > |m (|)| 2 X[i/i2,i/6]uhi/i2,-i/6](^i)- (52) 
On the other hand, by definition of mi, 

|mi(£i)| 2 = ( S mK 1 )) 2if E( if " 1 + 1( cos ^)) 2n 

n=0 V n / 

= (sin( 7 rei)) 2 ^'Hsm(7rei))^' £ " 1 + ") (cos^i))^ 

n=0 ^ ' 

Therefore, 

|mi(6)| 2 < min{l,C 2 |<i| 2 ^-^} (53) 

since |mi(£i)| 2 < 1. Let now 



g = 4vr(C 2 ) 1 /( 2 (^^)) ) g > = 2n l^ C2 Q | mo 



i 

Jl-l /«,_,• \ 2 \ 2 7 



3 2- J i +1 EJM")IH 



and r = 2q', 



with C*2 being defined in ( !21~]) . It is easy to check that g > r > g' > 0. By Propositions 14.21 
and S3] ([52]) and ([53]), we obtain (123"]) and 

In ( 123]) . the decay rate 7 is given by 7 = /<" — ZT' — |log 2 (C 2 ). Hence, by Lemma 143] in 
the special case that K' = and L > 10, we have 7 > L/4 + 1/2 > 3. This implies the 
feasibility condition ([6]). 
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